arXiv:1507.08066vl [astro-ph.GA] 29 Jul 2015 


Astronomy & Astrophysics manuscript no. GES'Submitted printer © ESO 2015 

July 30, 2015 


The Gaia-ESO Survey: Characterisation of the [a/Fe] 
sequences in the Miiky Way discs* 

G. Kordopatis^, R.F.G. Wyse^, G. Gilmore^, A. Recio-Blanco^, R de Laverny^, V. Hill^, V. Adibekyan^, U. Heiter®, 

I. Minchev^, B. Famaey^, T. Bensby®, S. Feltzing®, G. Guiglion^, A.J. Korn®, S. Mikolaitis^’®, M. Schultheis^, 

A. Vallenari^®, A. Bayo^^, G. Carraro^^, E. Flaccomio^®, E. Eranciosini^®, A. Hourihane®, P. Jofre®, S. E. Koposov®’^^, 
C. Eardo^®, J. Eewis®, K. Eind®, E. Magrini^®, E. Morbidelli^®, E. Pancino^^’^®, S. Randich^®, G.G. Sacco^®, 

C.C. Worley®, and S. Zaggia^® 

(Affiliations can be found after the references) 

Preprint online version: July 30, 2015 


ABSTRACT 

Context. High-resolution spectroscopic surveys of stars indicate that the Milky Way thin and thick discs follow different paths in the chemical 
space defined by [a/Fe] vs [Fe/H], suggesting possibly different formation mechanisms for each of these structures. 

Aims. We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined 
chemically as the high—a and low—a populations), and discuss their compatibility with discs defined by other means such as metallicity, kine¬ 
matics or positions. 

Methods. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chem¬ 
ical space for stars at different ranges of Galactocentric radii and heights from the Galactic mid-plane. The separation we find in velocity space 
allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without 
making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, adopting the 
magnesium abundance as a tracer of the a—elements, we characterise the spatial variation of the slopes of the [a/Fe] — [Fe/H] sequences for the 
thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. 

Results. We find that the thick disc, defined as the stars tracing the high—a sequence, extends up to super-solar metallicities ([Fe/H] « -1-0.2 dex) 
and the thin disc, defined as the stars tracing the low—a sequence, at least down to [Fe/H] « —0.8 dex, with hints pointing towards even lower 
values. Radial and vertical gradients in a—abundances are found for the thin disc, with mild spatial variations in its [a/Fe] — [Fe/H] paths, whereas 
for the thick disc we do not detect any such spatial variations, in agreement with results obtained recently from other high-resolution spectroscopic 
surveys. 

Conclusions. The small variations in the spatial [a/Fe] — [Fe/H] paths of the thin disc do not allow us to distinguish between formation models 
of this structure. On the other hand, the lack of radial gradients and [a/Fe] — [Fe/H] variations for the thick disc indicate that the mechanism 
responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the 
(present) thin disc started forming. 

Key words. Galaxy: abundances. Galaxy: disk. Galaxy: stellar content. Galaxy: evolution. Galaxy: dynamics and kinematics, stars: abundances 


1. Introduction 

The decomposition and separation of the Galactic disc into 
thin and thick counter-parts, first suggested by Gilmore & Reid 
(1983), has many consequences in terms of Galaxy formation 
and evolution. Indeed, the differences in chemistry and kinemat¬ 
ics between the stars of these two structures put constraints on 
both the role of mergers in the Galaxy in the last 8—12 Gyr 
(e.g.: Wyse 2001; Abadi et al. 2003; Villalobos & Helmi 2008) 
and on the importance of internal evolution mechanisms such 
as radial migration, dissipational cooling dependent on metal¬ 
licity or clumpy turbulent discs (e.g.: Sellwood & Binney 2002; 
Schonrich & Binney 2009; Minchev et al. 2013; Bournaud et al. 
2009). 

Among the literature, there is no unique way of selecting 
stars belonging to either of the discs. Often, kinematic criteria 
are used, but selections based on chemical composition or spatial 


* Based on observations collected with the FLAMES spectrograph 
at the VLT/UT2 telescope (Paranal Observatory, ESO, Chile), for the 
Gaia-ESO Large Public Survey, programme 188.B-3002. 


position are also employed, introducing potential biases in the 
distribution functions of the studied structures. Indeed, the first 
disc divisions, mainly based on star counts and density profiles, 
have found that the properties of the discs have an important 
overlap in kinematics and chemistry (Gilmore et al. 1989). 

The advent of multi-object spectroscopy and large spectro¬ 
scopic surveys such as SEGUE (Yanny et al. 2009) or RAVE 
(Steinmetz et al. 2006), confirmed these facts and measured 
the different kinematics relations of the discs (e.g. Nordstrom 
et al. 2004; Kordopatis et al. 2011, 2013c) and metallicities 
([Ee/H]thin ^ -0.1 dex, [Ee/H]thick ~ -0.5 dex, e.g: 
Bensby et al. 2007; Reddy & Lambert 2008; Ruchti et al. 2011; 
Kordopatis et al. 2013b). However, despite some specific at¬ 
tempts to characterise the wings of the distributions (e.g.: Norris 
et al. 1985; Chiba & Beers 2000; Carollo et al. 2010; Ruchti 
et al. 2011; Kordopatis et al. 2013b, 2015) the overlap in kine¬ 
matics and metallicity for the discs implies that the shape and 
extent of their distribution functions remain still quite uncertain. 
The best discriminant parameter seems to be the stellar age, the 
thick disc stars being consistently older than 9 Gyr, (e.g.; Wyse 
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& Gilmore 1988; Edvardsson et al. 1993; Bensby et al. 2005; 
Bergemann et al. 2014; Haywood et al. 2013). Yet, stellar age 
is a parameter difficult to obtain, requiring either asteroseismic 
data, or parallax measurements combined with high-resolution 
spectroscopy. 

In addition to age, the different formation histories of the two 
discs seems also strongly supported from the individual abun¬ 
dances measurements obtained for very local volumes of low 
mass FGK stars from very high-resolution spectra (R > 35 000, 
e.g.; Fuhrmann 1998, 2008, 2011; Reddy et al. 2006; Adibekyan 
et al. 2012; Bensby et al. 2014). In particular, the study of 
Bensby et al. (2014) used kinematic criteria to define the thin 
and thick discs, and found that each of these structures fol¬ 
lowed different sequences in the parameter space formed by the 
abundance ratio of a—elements over the iron, as a function of 
iron ([a/Fe] — [Fe/H]), with potentially a gap between these 
sequences. These results support previous conclusions that the 
thick disc stars were formed on a short time-scale, with rel¬ 
atively little pollution of the interstellar medium from super¬ 
novae of type la which are the main source of iron (e.g:. Wyse 
& Gilmore 1988; Freeman & Bland-Hawthorn 2002, and refer¬ 
ences therein). 

All of the previously cited high-resolution surveys were, 
however, limited in volume (typically only probing up to ~ 
100 pc from the Sun) and often with complex selection functions 
combining metallicity and kinematics. Further high-resolution 
spectroscopic investigation was therefore needed in order to as¬ 
sess whether this separation was a global feature in the discs or a 
result of an observational bias. The gap could indeed disappear 
or become a trough farther away from the Sun, in which case the 
properties of the stars with intermediate a—abundances would 
shed important light on the way the disc evolved chemically and 
dynamically. 

Recio-Blanco et al. (2014) used the first data-release (iDRl) 
of the Gaia-FSO Survey (Gilmore et al. 2012), analysing faint 
FGK stellar populations with high resolution spectra (R ^ 
20 000) from the FFAMFS/GIRAFFE spectrograph, to identify 
these two distinct sequences out to much larger volumes (up to 
distances of 5 kpc from the Sun) and made a first attempt to de¬ 
scribe their chemo-dynamical characteristics. This analysis was 
pushed further by Mikolaitis et al. (2014), who derived, for the 
same observed dataset, abundances of eight elements (Mg, Al, 
Si, Ca, Ti, Fe, Cr, Ni, Y). Mikolaitis et al. (2014) concluded, con¬ 
sistent with other analyses, that the gap between the sequences 
is the most significant when considering the [Mg/Fe] abundance 
ratio as a function of [Fe/H], and that the two discs show no 
significant intrinsic dispersion in their a—abundances at fixed 
[Fe/H]. This lack of scatter in the abundances along the metallic¬ 
ity for both thin and thick disc indicates first of all that the initial 
mass function (IMF) of the stars contributing to the [a/Fe] ratio 
at low metallicities and the fraction of binary stars at the ori¬ 
gin of supemovae type la, was constant over the formation time 
of the stars they pre-enrich, and in addition that the interstellar 
medium (ISM) should be very well mixed in order to have the 
same (small) dispersion at all metallicities. 

Using only data for Red Clump stars from the Apache 
Point Observatory Galactic Evolution Experiment (APOGEE), 
Nidever et al. (2014) measured the [a/Fe] — [Fe/H] trend for 
the thick disc and found that it had only very small variations 
with Galactocentric radius. Comparing with the results of a sim¬ 
ple Galactic chemical evolution model, the authors concluded 
that the small amplitude, less than 10 per cent, of the spatial 
variations is possibly indicative of an essentially constant star 
formation efficiency, and constant outflow rate balancing the 



Fig. 1. Aitoff projection of the Galactic coordinates of the Gaia-ESO 
iDR2 targets used in this study. 


Star-formation rate, during the early stages of the Milky Way 
while the stars that now form the thick disk were formed. Their 
favoured model has a star-formation efficiency corresponding to 
a gas consumption timescale of 2 Gyr, close to that inferred for 
clumpy discs observed at redshift 2 (e.g. Tacconi et al. 2013). 

We here use the FFAMFS/GIRAFFE Gaia-ESO second in¬ 
ternal data-release (iDR2) with its increased statistics (roughly 
a factor of two on the number of observed targets compared to 
iDRl), its enhanced analysis (use of a wider wavelength range in 
the spectral analysis due to an improvement of the parameter de¬ 
termination pipelines) and better quality (higher signal-to-noise 
ratio for the spectra, on average). We aim to characterise in a ro¬ 
bust statistical way the chemical paths followed by the thin and 
thick discs, in velocity and space, in order to provide a thorough 
insight on its possible origins. The understanding of the chemi¬ 
cal path origins will eventually further constrain disc formation 
theories, by indicating, for example, whether the thick disc was 
formed entirely in situ or by the contribution of extragalactic 
material (gas and/or stars). 

Section 2 describes the dataset and the quality cuts that have 
been used in this work. Section 3 identifies in a novel man¬ 
ner the presence of the two disc populations by showing the 
way a—abundances correlate with velocities. Section 4 identi¬ 
fies the two populations in terms of star-counts and investigates 
how their properties change at different Galactic regions. Finally, 
Sect. 6 summarises and concludes the present study. Note that in 
this paper we use the term metallicity to designate both the ratio 
of the global metal abundance to that of hydrogen ([M/H]) and 
the ratio of the true iron abundance to that of hydrogen ([Fe/H]). 
In general, the abundance pipeline for Gaia-ESO survey pro¬ 
vides results such that [M/H] « [Fe/H] (Recio-Blanco et al., 
in prep.), but we make explicit mention of one or the other pa¬ 
rameter when confusion could occur. 


2. Description of the data: Gaia-ESO iDR2 

We use throughout this paper the parameters derived from spec¬ 
tra obtained with the FFAMFS/GIRAFFE HR 10 and HR21 set¬ 
tings (see Fewis et al., in preparation, for a description of the 
data reduction procedures). The Galactic coordinates (f, b) of the 
observed targets are illustrated in Fig. 1. The source catalogue 
from which these targets are selected is the VISTA Hemisphere 
Survey (VHS, McMahon et al in prep.), with the primary selec¬ 
tion being: 

• 0.0 <{J -K)< 0.45 for 14.0 < J < 17.5 

• 0.4 < {J-K) < 0.70 for 12.5 < J < 15. 
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Extinction was taken into account by shifting the colour- 
boxes by 0.5 X E{B — V), where E{B — V) is the median 
Schlegel et al. (1998) extinction in that line-of-sight. If not 
enough targets were available within the colour cuts, additional 
targets were assigned by relaxing the red-edge of the colour-cut 
(see Gilmore et al. 2012, and Gilmore et al. in prep, for further 
details). 

We adopted the recommended parameter values from the 
Gaia-ESO analysis for the effective temperature (Tes), surface 
gravity (log g) and elemental abundances and the iron abundance 
we use is derived from individual Ee I lines. The determination 
of the atmospheric parameters for the Gaia-ESO survey is de¬ 
scribed in Recio-Blanco et al. Eor the sake of completeness we 
give here a brief description of how these were obtained and re¬ 
fer the interested reader to the original reference for further de¬ 
tails. We note that throughout the paper, the signal-to-noise ratio 
(S/N) refers to the value achieved in the HRIO setup (centred at 
548.8 nm). Eor a given star, the S/N in the HR21 setting (centred 
at 875.7nm) is approximately twice as large. 

2.1. Determination of the atmospheric parameters and 
abundances 

The stellar atmospheric parameters and abundances were de¬ 
termined by the members of the Gaia-ESO work package in 
charge of the GIRAEEE spectrum analysis for EGK-type stars. 
The individual spectra were analysed by three different proce¬ 
dures: MATISSE (Recio-Blanco et al. 2006), FERRE (Allende 
Prieto et al. 2006, and further developments) and Spectroscopy 
Made Easy (SME, Valenti & Piskunov 1996, and further devel¬ 
opments). 

Each of these procedures performed well with respect to 
APOGEE results and well-studied open and globular clusters 
(i.e. the comparisons resulted in low-amplitude dispersions and 
simple biases). The random errors on the parameter values were 
reduced by putting the results from each procedure on the same 
scale, then combining them. This was achieved by first choos¬ 
ing the results of a given procedure as the reference, and then 
multi-linear transformations as a function of Tes and log g were 
used to put the results from the two other procedures on the same 
scale as the reference one. The adopted transformations were of 
the form: 

= Op Oj- = cq Cl • TeffjP -f C 2 • log gp (1) 

where 0 is one of Tes, log (/, [M/H] or [a/Fe], the subscript r 
indicates the value from the reference procedure, and the sub¬ 
script p indicates the value from either of the other procedures 
(see Recio-Blanco et al., in prep., for more details). 

The resulting random errors for each method were fairly sim¬ 
ilar (~ 60K for Tgff, 0.14 for log;/, 0.06 for [M/H] and 0.05 
for [a/Fe]) and, as a consequence, the results for each star were 
combined by a simple average with the same weight given to 
each method. The averaged results were then calibrated with re¬ 
spect to 19 benchmarks stars (Jofre et al. 2014); this corrected for 
systematic biases of -1-83K for Tes and -1-0.23 dex for [M/H]. 

Once the final estimates of the values of the atmospheric 
parameters were obtained, the individual abundances, includ¬ 
ing the magnesium (Mg) and iron (Fe) abundances used in this 
work, were derived from three independent methods: an au¬ 
tomated spectral synthesis method (Mikolaitis et al. 2014), a 
Gauss-Newton method using a pre-computed synthetic spectra 
grid (Guiglion et al. 2014) and SME (Valenti & Piskunov 1996). 
In line with the Gaia-ESO consortium analysis, all the abun¬ 
dance analysis methods used MARCS (Gustafsson et al. 2008) 


model atmospheres. In addition, the input atomic and molecu¬ 
lar data were provided by the linelist group of Gaia-ESO, who 
collated the most recent and complete experimental and theoret¬ 
ical data sources. The comparison of the results of each method 
showed that a simple shift was enough to put them all on the 
same scale. This shift was independent of both the atmospheric 
parameters and abundances, and corrected for small offsets with 
respect to the Solar abundance (of order ±0.05 dex). Finally, for 
each star and each element, the results of the three methods were 
averaged with equal weight, and the elemental abundances rel¬ 
ative to the Sun (denoted [X/Fe] for element X) were obtained 
using the Grevesse et al. (2010) Solar elemental abundances. 

The final mean (median) errors, corresponding to the param¬ 
eter dispersion within the analysis nodes that have determined 
the atmospheric parameters and abundances are: 70 K (36 K) 
for Teff, 0.14 dex (0.08 dex) for \ogg, 0.08 dex (0.06 dex) for 
[Fe/H] and 0.06 dex (0.05 dex) for [Mg/Fe]. 

2.2. Determination of Gaiactocentric positions and 
veiocities 

The absolute magnitudes of the stars have been obtained by 
projecting the atmospheric parameters (Tgft, iogg, [M/H]) 
and VISTA (J — Kg) colour of the stars on the Yonsei-Yale 
isochrones (Demarque et al. 2004), as described in Kordopatis 
et al. (2011, 2013c), with the updates made in Recio-Blanco 
et al. (2014) concerning the treatment of the reddening. The 
derived distances have a mean uncertainty of 15 per cent. 
They were combined with the Galactic coordinates, line-of- 
sight velocities of the stars and the proper motions from the 
PPMXL catalogue (Roeser et al. 2010) in order to compute the 
Gaiactocentric positions and velocities in a cylindrical frame, 
with Vr,V 0 and Vz defined as positive with increasing R,(j) 
and Z (with the last towards the north Galactic pole). Typical 
errors of the proper motions are 4 — lOmasyr”^ (Roeser 
et al. 2010 ), resulting to transverse velocity errors at 2kpc of 
37 — 94kms“^. The adopted Solar motion with respect to the 
Local Standard of Rest (LSR) is that of Schonrich et al. (2010), 
namely U© = 11.1, V© = 12.24, W© = 7.25kms"!. Finally, 
the LSR is presumed to be on a circular orbit with azimuthal 
velocity Vc = 220kms“^, and we take the Sun to be lo¬ 
cated at (Ro,Zo) = (8, 0) kpc. The associated errors on the 
Gaiactocentric positions and velocities have been obtained by 
running 5000 Monte-Carlo realisations on the distances, proper 
motions and radial velocities and by measuring the dispersion of 
the derived parameters. 

The mean stellar velocities (after the quality selections, see 
below) as a function of Gaiactocentric position are illustrated 
in Fig. 2, with a binning of 0.4 kpc. The maps show consistent 
results with those already published in the literature. In partic¬ 
ular, we detect a mild radial gradient in the velocity compo¬ 
nent along the radial direction Vr, consistent with that found 
by Siebert et al. (2011b). This is visible at negative Z, going 
from ~ 2 kms“^ at i? ^ 6.5 kpc to ~ —5 kms“^ at 9 kpc (Vr, 
left panel of Fig. 2). We also detect a compression/rarefaction 
pattern in the vertical velocity component, seen as the difference 
in sign of < Vz > at Z = ±lkpc and R = 6, 9 kpc (right 
panel. Fig. 2), consistent with the patterns identified in SEGUE 
G-dwarf data by Widrow et al. (2012) and in RAVE data by 
Williams et al. (2013) and Kordopatis et al. (2013b). These mo¬ 
tions may be induced by external perturbations such as a satellite 
galaxy (e.g. Widrow et al. 2014) or internal perturbations such 
as spiral arms (e.g. Faure et al. 2014). 
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Fig. 2. Mean radial (left), azimuthal (middle) and vertical (right) Galactocentric velocities of the FLAMES/GIRAFFE stars in the Gaia-ESO 
internal Data-Release 2 (iDR2) as a function of their position in the Galaxy (in a cylindrical reference frame). For each panel, only the stars having 
errors less than 65kms“^ in that particular component of space velocity are selected, resulting to 3162, 3667, 3243 targets, for Vr, and Vz 
respectively. Mean values have been computed in bins of 0.4 kpc size, provided each position-bin contained at least 15 stars, resulting to an error 
on the mean velocity of 5 km per bin. 
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Fig. 3. Scatter plot of the iron abundances of the stars having each of 
the components of their Galactocentric space velocity known better 
than 65kms“'^ as a function of their measured Galactocentric radius 
and distance from the plane. The position of the Sun, at {Ro, Zo) = 
(8, 0) kpc, is indicated with a black ‘-I-’ symbol. The dashed box rep¬ 
resent the selection made for the “local” sample around the Sun (see 
Sect. 3 and Sect. 4.4). Typical error on positions is 8 per cent. 


2.3. Selection of subsamples based on quality cuts 

The Gaia-ESO Survey observes three classes of targets in the 
Galaxy (field stars, members of open clusters and calibration 
standards, see Gilmore et al. 2012). Our analysis in this paper 
focuses on the Milky Way discs and therefore we kept only the 
field stars, i.e. those stars that are neither members of open clus¬ 
ters nor part of the calibration catalogue. 

Furthermore, in order to have a dataset with reliable param¬ 
eters, we excluded stars with spectra having a signal-to-noise 
ratio (S/N) less than 10 in the HRIO setup (higher S/N cuts have 
also been tested, without changing the overall conclusions of this 
paper) or having the highest uncertainties in the values of the de¬ 
rived parameters (i.e. being at the tails of the error distributions). 
As a compromise between accuracy and sample size, we set the 
cuts for the errors in effective temperature at crTeft ^ 250 K 


(4 per cent of the stars), in surface gravity at uiogg > 0.5 (3 per 
cent of the stars), in iron abundance at (T[Fe/H] > 0-3 dex (0.5 per 
cent of the stars), and in line-of-sight distance at ao/D > 0.5 
(4 per cent of the stars). 

In addition, for each velocity component studied, we ex¬ 
cluded those stars having errors in the given velocity greater 
than 65kms“^. The errors in three-dimensional space veloc¬ 
ity are dominated by the errors in the tangential velocities, so 
these rejection criteria lead to selection functions for each ve¬ 
locity component that are dependent on the stellar distances and 
the particular lines-of-sight (see Fig. 3). Indeed, stars close to 
the Galactic plane tend to have large errors on Vz, and stars 
towards the Galactic poles have poorly constrained Vr. The re¬ 
moval of stars that had a derived error in a given velocity com¬ 
ponent that was greater than 65kms“^, left us with 3162, 3667 
and 3243 stars, for Vr, V<^ and Vz, respectively. The spatial dis¬ 
tribution of the 2797 stars that have errors in all three of their 
velocity components below 65kms“^ is illustrated on Fig. 3, 
colour coded according to iron abundance. One can see from 
this plot that the metal-rich population ([Fe/H] > —0.2 dex) is 
mostly found close to the plane, below 1 kpc, as expected for a 
thin-disc dominated population. 

Figure 4 shows the scatter plot of the azimuthal velocities 
of the 3667 stars that passed the error cut in this velocity com¬ 
ponent, in the chemical space defined by magnesium and iron. 
Following Mikolaitis et al. (2014), magnesium is used as the 
a—element tracer in the remainder of this paper. The gap sug¬ 
gested from the analyses of Adibekyan et al. (2012, using the 
HARPS spectrograph), Bensby et al. (2014, using mainly the 
MIKE spectrograph) and Recio-Blanco et al. (2014, using the 
high S/N sample of Gaia-ESO iDRl) is seen in Fig. 4 as a trough, 
due to the low S/N threshold (S/N> 10) we imposed for the 
analysis at this stage (median S/N= 25; see Sect. 4.2 for an 
analysis of mock data illustrating how a real gap can get pop¬ 
ulated by noisy data and Sect. 4.4 for an analysis using higher 
S/N cuts on the Gaia-ESO iDR2 data). Indeed, the two popula¬ 
tions - one low—a ([Mg/Fe] < 0.2 dex) on orbits with higher 
azimuthal velocities (mean V^ ~ 210 kms“^), and one high—a 
([Mg/Fe] > 0.1 dex) on orbits with lower azimuthal velocities 
(V^ ^ 185kms“^) - can be identified with the colour-code. 
These two populations are more clearly seen in Fig. 5, repre¬ 
senting V^ as a function of [Fe/H], colour-coded according to 
the [Mg/Fe] abundance. The fact that the high—a population has 
a larger velocity dispersion is also evident. 
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Fig. 4. [Mg/Fe] as a function of [Fe/H] for the Gaia-ESO iDR2 stars 
with S/N> 10. The mean error bar is represented in grey at the bot¬ 
tom left comer. The range of the plot has been truncated to —1.5 < 
[Fe/H] < 0.55 dex. The contour lines are drawn for 33 per cent, 66 per 
cent, 90 per cent and 99 per cent of the sample. 
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Fig. 5. Azimuthal velocity as a function of [Fe/H[, colour-coded ac¬ 
cording to the [Mg/Fe[ abundance (the velocity and metallicity ranges 
have been tmncated for visualisation purposes). Dashed horizontal lines 
are plotted at = 210 kmand V<f, = 185 km. Typical errors 
on and [Fe/H] are represented at the bottom left comer, in grey. 

Given their different global kinematic properties, throughout 
this paper the high—a and low—a populations will be referred 
to as the thick and thin discs, respectively. In the next section, we 
explore the metallicity ranges up to which these two populations 
extend, by investigating their kinematics, but without decompos¬ 
ing the stellar populations. 

3. The local disc dichotomy in terms of 
Galactocentric velocities 

In this section we aim to investigate, for the stars in the greater 
Solar neighbourhood (7.5 < R < 8.5 kpc , |^| < Ikpc) - the 
Solar suburbs - the metallicity range over which the disc(s) can 
be decomposed into at least two different populations, defined by 


their correlation between kinematics and chemical abundances. 
This we achieve with no a priori assumption about the shapes 
of their metallicity distribution functions (DF) or their kinematic 
DFs. Kinematics provide the means to discriminate among mod¬ 
els for the formation and evolution of the discs, and the way in 
which metallicity and chemical abundances correlate with kine¬ 
matics carries important and fundamental information on the 
nature of the stars in our sample (e.g. Eggen, Lynden-Bell, & 
Sandage 1962; Freeman & Bland-Hawthorn 2002; Rix & Bovy 
2013). For the remainder of this paper, we dehne a stellar popu¬ 
lation by demanding that its member stars show the same corre¬ 
lations between chemistry and kinematics over the whole range 
of metallicity that the population spans. The approach presented 
below is complementary to, and independent of, that of Sect. 4, 
where we ht the a—abundances in different metallicity bins in 
order to trace the evolution of the trends for both discs, dehned 
purely by chemistry. 

3.1. The greater Solar neighbourhood sample 

The disc stars in the greater Solar neighbourhood are selected 
by requiring distances close to the Galactic plane {\Z\ < 1 kpc), 
Galactocentric radii between 7.5 and 8.5 kpc {i.e. distances 
from the Sun closer than 1.12 kpc) and azimuthal velocities 
greater than 0 km s“^. These selections ensure that; 

(i) The radial and vertical gradients in the discs have 
minimal effects on the interpretation of the results. 

(ii) The sample under consideration is as homogeneous as 
possible towards all the lines-of-sight (see Fig. 3). Excluding 
stars farther than 1 kpc from the plane should not particularly 
bias the kinematics of the selected populations (other than the 
bias coming from possible vertical gradients in V^), provided 
that the populations are quasi-isothermal in vertical velocity 
dispersion, which appears to be a reasonable approximation 
(e.g.: Elynn et al. 2006; Binney 2012; Binney et al. 2014). 

(Hi) The contamination by inner halo stars is reduced, in 
order to focus on the thin and thick disc chemical dichotomy. 
The adopted azimuthal velocity threshold of Okms“^ allows 
the retention of metal-weak thick disc stars, given an expected 
azimuthal velocity distribution for the thick disc that extends 
down to at least 60 — 80kms“^ (i.e. ^ 0kms“^ with a 
65kms“^ error, see Kordopatis et al. 2013b,c). This threshold 
simultaneously excludes counter-rotating stars likely to be 
members of the halo (in any case, the contribution from the halo 
should be small, given the metallicity range and distances above 
the plane investigated; see for example, Carollo et al. 2010). 
Note that in this work we treat the metal weak thick disc as part 
of the canonical thick disc; whether or not the chemo-dynamical 
characteristics of the metal-weak thick disc are distinct will be 
investigated in a future paper. 


3.2. Description of the method 

We assess the presence of two distinct populations in chemical 
space by measuring the mean velocities (Vi) and their disper¬ 
sions (uvi) for different [Fe/H] bins, in addition to the correla- 
tion(s) between the means of each of the velocity components 
and the stellar a—enhancements (noted simply i9Vi/£)[Mg/Fe]). 
While the evolution of the mean velocities as a function of metal- 
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licity has long been used to identify the transition between the 
thin and thick discs (e.g.: Gilmore et al. 1989; Freeman & Bland- 
Hawthorn 2002; Reddy et al. 2006), here the novel approach is 
to see whether there is a difference in the mean velocities at a 
given [Fe/H], as a function of [Mg/Fe]. Non-zero values for 
9Vi/i9[Mg/Fe] can arise due to the following: 

1. The estimated values of the atmospheric parameters, dis¬ 
tances and velocities have correlated errors and therefore bi¬ 
ases. 

2. At least one population has a velocity gradient and/or chem¬ 
ical gradient as a function of Galactocentric radius and/or 
vertical distance, resulting in different kinematic properties 
in a given metallicity bin. 

3. Only one population has an intrinsic correlation between 
chemistry and kinematics, as has previously been proposed 
for the thick disc in previous studies (c/V^/c/pe/H] « 
50kms“^ dex“^, e.g: Kordopatis et al. 2011, 2013b). 

4. There are at least two different populations with distinct 
kinematics and chemistry. 

5. A combination of some or all of the above. 

Possibility (1) is not likely. Indeed, as noted earlier 
(Sect. 2.1), the calibration using benchmark stars removed the 
known biases in the stellar parameters, while Recio-Blanco et al. 
(in prep.) found no significant degeneracies between the param¬ 
eter determinations that could have led to artificial correlations 
between the parameters. Nevertheless, hidden biases may still 
exist. For example, a bias in metallicity, effective temperature 
or surface gravity could affect the final distance estimation and 
thus the final Galactocentric space velocity. Biases of less than 
0.1 dex in [M/H] or less than 100 K in Tea (the maximum pos¬ 
sible values for this analysis), could change the distances up to 
10 per cent (see Schultheis et al. 2015, their Table 1, for further 
details). At 1 kpc, this would translate to changes in velocities of 
~ 5 km (assuming proper motions of 10 masyr“^), which is 
too small to have a significant effect on our analysis. 

Concerning possibility (2), the amplitude of chemical gra¬ 
dients in the thin disc has been thoroughly discussed in the 
literature, with evidence from a variety of stellar tracers for 
mild radial and vertical metallicity gradients, of the order of 
9[Fe/H]/c)R « —0.04 to —0.07 dexkpc”^ and (9[Fe/H]/c)Z « 
—0.05 to —O.ldexkpc”^ (e.g. Gazzano et al. 2013; Boeche 
et al. 2013; Recio-Blanco et al. 2014; Anders et al. 2014). 
The thick disc shows no evidence for a radial metallicity gra¬ 
dient (e.g.: Cheng et al. 2012; Bergemann et al. 2014; Boeche 
et al. 2014; Mikolaitis et al. 2014) and at most a mild verti¬ 
cal metallicity gradient (Ruchti et al. 2010; Kordopatis et al. 
2011, 2013c; Mikolaitis et al. 2014). The [a/Fe] gradients in 
both discs are also small (< ± 0.02 dexkpc“^), with the actual 
amplitude and even the sign being still a matter of debate (see 
Sect. 5 of Mikolaitis et al. 2014). As far as the kinematic trends 
are concerned, neither the radial gradient in Vr (OVr/OR « 
3 kmkpc” , Siebert et al. 2011b) nor the “compression- 
rarefaction” patterns reported in Vz in both R and Z (Williams 
et al. 2013; Kordopatis et al. 2013a), are strong enough to con¬ 
tribute significantly within the volume of the greater Solar neigh¬ 
bourhood. A similar conclusion holds for V^, given the stud¬ 
ies of Lee et al. (2011) and Binney et al. (2014) based, respec¬ 
tively, on the SDSS-DR 8 (Aihara et al. 2011) and RAVE-DR4 
(Kordopatis et al. 2013a) catalogues for stars located roughly in 
the same volume as the present study. 

The difference in the chemo-dynamical signature between 
that due to a single population having an intrinsic correlation be¬ 


tween metallicity and kinematics, possibility (3) above, or that 
due to the superposition of several populations each with distinct 
correlations between metallicity and kinematics (possibility (4)) 
is more subtle. It can however be seen in the variation of the gra¬ 
dient of kinematics with chemistry, such as 9Vi/c)[Mg/Fe] with 
metallicity. Indeed, a single population, as defined at the begin¬ 
ning of this section, should exhibit no trend in 9Vi/c)[Mg/Fe] as 
a function of metallicity. A combination of more than one pop¬ 
ulation would produce variations of this quantity if the relative 
proportions of the populations changes with [Fe/Hj. Our inves¬ 
tigation of the extent of the superposition of the stellar popula¬ 
tions using this approach therefore differs from those previously 
published in the literature, since here we make no assumptions 
about the number of these populations that may be present, nor 
the shape of their distribution functions (see for example Bovy 
et al. 2012 b, for a decomposition of the disc into multiple mono¬ 
abundance populations). 

Figures 6 , 7, 8 and 9 show the results obtained for the az¬ 
imuthal (V^), radial (Vr) and vertical (Vz) velocity components 
of the stars. Each of these figures represent the mean number of 
stars per iron abundance bin^ (panels A), their mean radial dis¬ 
tance from the Galactic centre (panels B), their distance from 
the plane (panels C), their mean velocity and dispersion (panels 
D and E), and finally the correlation between the mean veloc¬ 
ity and the a—abundances (panels E and Eig. 9). They have been 
obtained with 100 Monte-Carlo realisations on positions, veloc¬ 
ities and iron abundance, and by keeping the iron bin boundaries 
at fixed values. Stars can therefore move between bins from one 
realisation to the next, as parameter values change with each re¬ 
alisation. 

After each Monte-Carlo realisation, the mean velocities and 
Icr velocity dispersions of the stars within each bin were mea¬ 
sured, provided that bin contained at least 15 stars. As far as the 
()Vi/c)[Mg/Fe] values are concerned, they have been obtained 
by fitting a line in the Vi vs [Mg/Fe] parameter space and by 
measuring the slope of that line. The final values inside each 
bin were then obtained by computing the average over the re¬ 
alisations. The associated error bars were obtained by dividing 
the standard deviation of the means by the square root of the 
mean number of stars inside each bin. Finally, the corrected ve¬ 
locity dispersions (green lines in Figs. 6 , 7 and 8 ) were obtained 
by removing quadratically two times the mean velocity error of 
the stars inside each bin, the factor of two arising since the er¬ 
rors are taken into account twice: first, intrinsic to the measure¬ 
ment, and then again during the Monte-Carlo realisations (see 
also Kordopatis et al. 2013b). We therefore used the following 
relation to correct the velocity dispersions for observational un¬ 
certainties: 


avi = y <^'v- ~ 2 X mean{errorViY (2) 

where avi represents the corrected velocity dispersion of 
Vr,V 0 or Vz, and ay. is the corresponding measured veloc¬ 
ity dispersion (averaged over the Monte-Carlo realisations). We 
refer the reader to Guiglion et al. (submitted) for a thorough anal¬ 
ysis of the velocity dispersion of the disc stars as a function of 
the chemical composition. 


’ Figures 6 , 7, 8 each contain different number of targets, accord¬ 
ing to the selection criteria performed based on the velocity errors. The 
distributions (panels A) are not corrected for the selection function. 
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Fig. 6. Mean number of stars (A), Galactocentric radii (B), distances from the Galactic plane (C), azimuthal velocities (D), azimuthal velocity 
dispersions (E) as measured from the data (in hlack) and corrected hy the error measurements (in green), at different metallicity bins, for stars 
located within the greater Solar neighbourhood. Panel F shows the correlations between the mean azimuthal velocity and the a—abundance ratios 
for different metallicity bins. The error bars represent the 3cr Poisson noise averaged over the Monte-Carlo realisations. The bins are 0.07 dex 
wide, chosen as a compromise between the number of stars inside each bin and an adequate coverage of the metallicity range. 


3.3. Azimuthal velocity component 

The identification of the thin and thick disc through the analysis 
of the stellar kinematics is most cleanly achieved using the az¬ 
imuthal velocities of the stars, because of the asymmetric drift in 
the mean kinematics. Indeed, from the trends of the mean ve¬ 
locity (Fig. 6 , panel D), one can identify the thin disc component, 
with ^ 205 km s“^, dominating the star counts for metallic- 
ities above —0.4 dex. The thick disc, virtually unpolluted by the 
halo (due to the imposed cuts: jZj < 1 kpc and 14^ > 0kms“^), 
is distinguishable below metallicities of —0.5 dex, where there 
is a decrease in the mean azimuthal velocities. The velocity dis¬ 
persion (Fig. 6 , panel E) at metallicities above —0.2 dex has 
the value ~ 25kms“^. Below —0.5 dex, this increases to 
(70 ~ 45kms“^. These values agree with the estimates of the 
velocity dispersions of the old thin and thick discs as measured 
far from the Solar neighbourhood (Kordopatis et al. 2011), those 
measured for stars in the RAVE survey within roughly the same 
volume (Binney et al. 2014), and other values available from the 
literature (e.g.: Soubiran et al. 2003; Lee et al. 2011). 

As far as the correlations between and [a/Fe] are 
concerned, (Fig. 6 , panel F and Fig. 9, top plot), the be¬ 


haviour of c)V 0 / 9 [Mg/Fe] with increasing metallicity shows 
three regimes: 

- a value close to zero at super-solar metallicities, 

- a non-zero value of c)V 0 /£)[Mg/Fe] « — 120 kms“^ dex“^ 
below —0.5 dex, compatible with little variation within the 
errors, 

- an intermediate regime between metallicities 
[—0.5,-FO.2] dex, where the amplitude of d\d\M.g/¥e] 
monotonically increases with metallicity. We interpret this 
variation as indicative of a varying mixture of populations. 

We note that these trends are not depending on the adopted 
binning in metallicity, since they are also recovered when adopt¬ 
ing, for example, overlapping bins 0.07 dex wide, spaced by 
0.05 dex. Following from the discussion at the beginning of 
Sect. 3, these three regimes suggest that the relative propor¬ 
tions of the different stellar populations present in our sample 
do not vary significantly above + 0.2 dex, since we find no trend 
of dVrf,/ i9[Mg/Fe] with metallicity in this range. Similar conclu¬ 
sions can be made below —0.5 dex, although the larger error bars 
in this regime could hide weak trends. In terms of separate thin 


7 

































G. Kordopatis et al.: Characterisation of the [a/Fe] sequences in the Milky Way discs 



-1.0 -0.5 0.0 0.5 

[Fe/H] (dex) 


Fig. 7. Same as Fig. 6 but for the radial velocity component, Vr. 


and thick discs^, the trends we find imply that the metallicity dis¬ 
tribution of the thick disc extends at least up to Solar values, and 
that of the thin disc goes down to at least —0.5 dex. This over¬ 
lap in metallicity was already suggested from the a—abundance 
trends derived from very-high resolution studies of stars in the 
immediate Solar neighbourhood (d < 100 pc, e.g. Reddy et al. 
2006; Adibekyan et al. 2013; Bensby et al. 2014; Nissen 2015), 
and by Nidever et al. (2014) for stars in a wide range of Galactic 
radii (0 < i? < 15 kpc). 

We can further assess the presence of thin disc stars be¬ 
low [Fe/H] = —0.5 dex by making a rough separation of the 
stars within a given [Fe/H] bin into either high—a ([Mg/Fe] > 
0.3 dex) or low—a ([Mg/Fe] < 0.2 dex, see also Sect. 5.1) can¬ 
didates. This results in the low—a population having a median 

~ 212kms“^, whereas the median of the high—a pop¬ 
ulation is ~ 165 km s“^. The low—a population has a lag com¬ 
pared to the LSR that is typical of the thin disc, and the lag 
for the high—a population is consistent with estimates for the 
canonical and metal weak thick disc, provided there is a corre¬ 
lation between and [Fe/H] of 50 km dex“^, as found by 
Kordopatis et al. (2013a). We therefore conclude that the thin 
disc metallicity distribution extends down to at least [Fe/H] = 

^ Note that the metal-weak thick disc is not treated separately from 
the canonical thick disc, and that its chemo-dynamical properties will 
be thoroughly investigated in a future paper. 


—0.8 dex. Below that metallicity, the number of stars is too small 
to allow us to derive a reliable estimate of the mean velocity. The 
relative proportion of thin disc stars at low metallicities will be 
discussed further in Sect. 4. 

A natural outcome of the existence of thin disc stars below 
[Fe/H] < —0.5 dex is that the non-zero 9V0/i9[Mg/Fe] value 
that we measure is at least partly due to the mixture of thin and 
thick disc. However, within the errors. Fig. 9 shows a flat be¬ 
haviour below [Fe/H] < —0.5 dex, suggesting a constant frac¬ 
tion of thick/thin disc. This result will also be found when sepa¬ 
rating the stars according to their a—abundances. Indeed, while 
we cannot eliminate the possibility that there exists an intrinsic 
correlation between and [a/Fe] for the thick disc stars, as 
suggested by Recio-Blanco et al. (2014), we note that the dif¬ 
ferences that we hnd in the median values for [Mg/Fe] ~ 
0.35dex and [Mg/Fe] ~ 0.1 dex (165kmand212kms“^, 
respectively), could be sufficient to explain the measured value 
of i9V^/c)[Mg/Fe] « —120kms“^ dex“^, provided a hxed 
fraction of thick disc stars over thin disc stars of approximately 
60 per cent. This ratio is in good agreement with the local thin 
disc / thick disc ratio of 3:2 for the metallicity bin [—0.8, —0.6] 
derived by Wyse & Gilmore (1995) through the combination of 
a volume-complete local sample {d < 30 pc) and an in situ sam¬ 
ple at Z ~ 1 kpc (after scaling with appropriate density laws 
and kinematics). 
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Fig. 8. Same as Fig. 6 but for the vertical velocity component, Vz- 


3.4. Radial and vertical velocity components 

Due to the fact that both the radial and vertical velocity distribu¬ 
tion functions are centred at zero for both discs^, the thin/thick 
disc dichotomy is harder to identify with these velocity compo¬ 
nents. The investigation of the mean velocities as a function of 
metallicity confirms indeed this fact (see panel D of Fig. 7 and 
Fig. 8), where we find mean velocities centred around 0 km s“^. 

The velocity dispersions (panels E) of the most metal- 
poor stars are found to be {api,az) ~ (50,40) km s“^, typi¬ 
cal of thick disc values, while the most metal-rich stars have 
dispersions more typical of the thin disc, being {an, az) ~ 
(30, 20) kms“^ (e.g.: Soubiran et al. 2003; Lee et al. 2011; Bovy 
et al. 2012a; Sharma et al. 2014). In particular, the abrupt change 
in the slope daz /9[Fe/H] for metallicities > —0.4 dex suggests 
that this is the metallicity threshold above which the thin disc be¬ 
comes the dominant population. This point will be conhrmed in 
Sect. 4.4, where the discs are separated in chemical space only. 

As far as the i9VR_z/9[Mg/Fe] trends are concerned, we 
hnd them being consistent with zero at least within the metal¬ 
licity range [—0.8, 4-0.2] dex. The results using showed that 
the relative proportion of the populations is changing across 


^ The compression/rarefaction patterns in the kinematics of the discs 
lead to only small deviations from zero means, see Siebert et al. (2011a); 
Williams et al. (2013); Kordopatis et al. (2013b). 


this metallicity range. Therefore, the invariant, null value for 
5VR_z/i9[Mg/Fe] implies that all populations present in our 
sample share common radial and vertical velocity distributions, 
at least in the mean. 

An investigation of the detailed trends of the velocity disper¬ 
sions is outside the scope of this paper, but we did investigate the 
specific case of the behaviour of daz /9[Mg/Fe] as a function of 
metallicity, as this may be compared with the independent analy¬ 
sis of Guiglion et al. (2015, submitted). We find a non-zero value 
for the metal-poor stars (though still marginally compatible with 
zero at the 3cr level). This value, suggests that in a given metal¬ 
licity range the thick disc stars have a larger velocity dispersion 
than the thin disc ones {i.e. the low—a ones), as expected (see 
Guiglion et al. for details). 

3.5. Discussion 

Our analysis of the mean azimuthal velocities of the stars shows 
that for metallicities between —0.8 and 4-0.2 dex, there is a mix¬ 
ture of at least two populations in our sample. Further, their rela¬ 
tive proportions vary with metallicity. This is a well-established 
result for the thin and thick disks defined by kinematics (e.g. 
Carney et al. 1989; Wyse & Gilmore 1995; Bensby et al. 2005, 
2007); the new result here refers to the discs when dehned by 
distinct chemical abundance patterns. 
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Fig. 9. Same as panels F of Figs. 6 to 8, hut put side hy side to facilitate 
comparison of the trends. The plots have been smoothed with a boxcar 
average taking into account the closest neighbour. 


Our results show that the thin disc, dehned as the 
low—a population, extends in metallicity at least as low as 
—O.Sdex, whereas the thick disc, dehned as the high—a pop¬ 
ulation, extends at least up to Solar metallicities. Outside the 
[—0.8, -I-0.2] dex range, the large error bars that we derive in our 
analysis cannot allow to draw any robust conclusions. 

The metallicity ranges over which we have found the discs to 
extend have to be considered in light of the results of Haywood 
et al. (2013) and Snaith et al. (2014). These authors concluded 
that chemical evolution was traced by the thick disk up to high 
metallicities, then after dilution by metal-poor gas, proceeded 
following the thin disk sequence from low metallicity to high 
metallicity. A thick disc extending up to super-Solar metallici¬ 
ties plus a thin disc extending down to —0.8 dex, challenge this 
view since such a large metallicity overlap - over 1 dex - would 
require accretion of a signihcant mass of metal-poor gas to lower 
the ISM’s metallicity before the onset of star formation in the 
thin disc. We note, however, that we have not corrected for the 
selection function in our analysis. Further, we have very limited 
information concerning stellar ages. 

Future data from Gaia will provide better age estimates 
from improved distances and absolute magnitudes, in addition 
to proper motions that are ten times more accurate than those we 
adopted here, which will allow a more refined analysis. 


4. Identification and characterisation of the 
[Mg/Fe] sequences in terms of star-counts 

In the previous section we identified and conhrmed the exis¬ 
tence of two populations in our Solar suburb sample, for metal¬ 
licities between [—0.8, -1-0.2] dex. In that analysis we made no 
assumptions about the underlying metallicity distribution func¬ 
tions. In this section, we extend the analysis to an investiga¬ 
tion of the a—abundance sequences across the entire metallic¬ 
ity range of the data. We aim to characterise the possible dou¬ 
ble structure in each metallicity bin, and obtain estimates of the 
mean a—enhancement value, and the dispersion, of each com¬ 
ponent, for different regions of the Galaxy (i.e. both within and 
beyond the extended Solar neighbourhood). As in the previous 
sections, we use the magnesium abundance as a tracer of the 
a—abundance of the stars, and we make the reasonable assump¬ 
tion that at a fixed [Fe/H] the Gaia-ESO survey is not biased 
with respect to the value of [Mg/Fe]. 

4. 1. Model fitting of the data 

The model-htting is done using a truncated-Newton method 
(Dembo & Steihaug 1983) on the un-binned [Mg/Fe] data of 
a given metallicity range, by maximising the log-likelihood of 
the following relation; 


N 2 

cri, Pi,/X2, cr2j P 2 ) = HE 


Pi 


e 


([a/Fe]t-Mi)^ 


(3) 

where N is the total number of stars considered in a given metal¬ 
licity range, pi, p. 2 j are the mean [a/Fe] values of the thin and 
thick disc, cri, (T 2 are the dispersions on the [a/Fe] abundances. 
Pi and p 2 are the relative proportions of each population, and 
[a/Fe]t is the [a/Fe] abundance ratio of a given star measured 
from its spectrum. The fit is performed for every Monte-Carlo 
realisation separately. 

Motivated by the results of Adibekyan et al. (2013); Bensby 
et al. (2014); Recio-Blanco et al. (2014), we assume that the 
mean [a/Fe] abundance of each population with [Fe/H] > 
—0.9 dex is a function of metallicity, and that the thin disc’s 
mean [a /Fe] at any metallicity is always lower than (or equal to) 
the mean [a/Fe] abundance of the thick disc. We note that the 
thick disc exhibits an [a/Fe] plateau at lower metallicities (e.g. 
Ruchti et al. 2010), this plateau ending with the so-called knee 
of the [a/Fe] — [Fe/H] distribution. The exact position of the 
knee is a matter of debate; its location depends on how quickly 
star formation and chemical enrichment happened in the thick 
disc. The quality of our data at the lowest metallicities does not, 
however, allow us to investigate this issue. 

To perform our htting, the initial assumed values for the 
mean a—enhancements of the discs are: 


/^thick — —0.30 X [Fe/H] + 0.2 (4) 

Pthin = -0.10 X [Fe/H]. (5) 

We allow the fitting procedure to vary these means freely by 
up to ±0.15 dex. As far as the dispersions of the distributions 
are concerned, we allow them to vary in the interval between 
[0.05,0.15] dex, and start with an initial value of 0.08 dex, larger 
than the typical errors in [a/Fe] (see Sect. 2.1). 

This set of assumptions allows the value of [Mg/Fe] at the 
low metallicity end of the thick disc distribution to be between 
0.35 and 0.65 dex (consistent with the results of Mikolaitis et al. 
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2014, obtained with the Gaia-ESO DRl dataset). Furthermore, 
the priors also ensure that the thin disc sequence has enhance¬ 
ments that are close to Solar values at Solar metallicities and 
that the [Mg/Fe] prior for the thick disc is always higher than 
that for the thin disc (the starting relations are equal at ten times 
the solar iron abundance). 

Note that our procedure always fits two components to the 
data, even in regions where only one sequence could exist. Since 
a model with two components will always fit the data better, we 
investigate the significance of this fit by deriving the p—value 
of the log-likelihood ratio between the null hypothesis of having 
only one thick/thin disc and the alternative model of having both 
components. We therefore define the following test statistic of 
three degrees of freedom: 

£> = 2 • (log £2 - log£i) (6) 

where Li is the likelihood associated to the model with two com¬ 
ponents (5 degrees of freedom, see Eq. 3) and £i is the likeli¬ 
hood associated to the null model with only one component (2 
degrees of freedom). We set that for a given Monte-Carlo realisa¬ 
tion the two components model is the more likely if > 7.82, 
which corresponds to a p—value < 0.05. We then derive the 
mean fitted populations over all the Monte-Carlo realisations, 
and conclude that two populations are needed only if the mean 
is greater than 1.5 {i.e. the two component model is preferred for 
more than half of the realisations with 0.95 confidence). 

This procedure has been applied to six subsamples, selected 
by location within the Galaxy: three ranges in Galactocentric ra¬ 
dius (inner Galactic ring: 6.5 — 7.5 kpc. Solar suburbs: 7.5 — 
8.5 kpc, outer Galactic ring: 8.5 — 10 kpc), and two distances 
from the plane (\Z\ < Ikpc and 1 < \Z\ < 2kpc). The re¬ 
sults, obtained after averaging 100 Monte-Carlo realisations on 
distances, metallicities and [Mg/Fe] abundances, are discussed 
in the following sections. 

4.2. Tests on mock data 

We tested our fitting procedure on mock data to verify that it is 
able to find the “correct” a— enhancements of two populations 
with chemical trends similar to what is expected in the Milky 
Way. For that purpose, we took the [Fe/H] values of the stars in 
our sample, and assigned for the range —0.8 < [Fe/H] < -1-0.3 
the trends seen on the top left plot of Fig. 10: the position of 
the thick disc knee is located at —0.8 dex, the thin disc’s lowest 
metallicity is —0.6, and the two sequences merge at [Fe/H] > 
0.3. These values are selected in order to roughly represent the 
actual data, while still differing from the priors of Eqs.4 and 5; 
this should allow us to verify that our procedure converges to¬ 
wards the correct results. The ‘stars’ were randomly assigned to 
either the high— a or low— a population. Then, for each ‘star’ 
we assigned the uncertainties of the true Gaia-ESO iDR2 cata¬ 
logue, to reproduce realistic error bars. The blurring due to the 
uncertainties on the iron and a— abundances can be seen in the 
top-right panel. One can notice that the underlying gap in the 
chemical paths is now barely visible and that Fig. 10 is qual¬ 
itatively similar to Fig. 4 where the actual Gaia-ESO data are 
plotted. 

We applied our procedure, starting from equations (4) and 
(5), to 100 Monte-Carlo realisations of noisy mock data and de¬ 
rived the average trends plotted in red crosses and green square 
symbols on the bottom left panel of Fig. 10. One can see that the 
“true” trends are nicely recovered, up to Solar metallicities, vali¬ 
dating our procedure. In particular, the procedure performs well 


for a gap-width varying with metallicity, even when the sepa¬ 
ration is very small. Our procedure also correctly finds that at 
the low metallicities a single high—a population is enough to 
fit the data. Indeed, despite fitting two populations even at the 
lowest metallicities (green squares below [Fe/H] = —0.6 in 
the lower left panel of Fig. 10), the relative weight associated 
to these points, combined with the test statistic of Eq. 6, indicate 
that these points are not likely to represent correctly the data (the 
red line in the bottom right plot of Fig. 10 is consistently below 
1.5). 

The difficulty of fitting the data towards the highest metal¬ 
licities, where the the separation between the two sequences is 
small, can be seen on the noisiness of the average log-likelihood 
ratio, comparing the two-component model with a model having 
only the thin disc (green line). We note, however, that the two 
component model is always the preferred fit, even at this regime, 
as it should be, based on the input mock data. 

In the following sub-sections, we therefore apply our proce¬ 
dure with confidence to the Gaia-ESO iDR2 data. 


4.3. Robustness to signal-to-noise cuts 

Figure 11 shows the trends of the high—a and low—a popu¬ 
lations at different Galactocentric radii and distances from the 
plane, colour-coded according to different cuts in S/N of the 
actual Gaia-ESO spectra. One can see that our results are ro¬ 
bust to these cuts, with the overall trends being similar, inde¬ 
pendent of the S/N threshold. Lowering the threshold increases 
the number of analysed stars and populates regions of chemi¬ 
cal space which statistically have fewer stars, i.e. at lower and 
higher metallicities. A lower S/N threshold therefore achieves 
a better description of the overall behaviour of the trends, but 
with larger uncertainties on individual estimates, particularly at 
the extremes of the metallicities. In the remainder of this analy¬ 
sis, when not stated otherwise, we adopt S/N> 25 as providing a 
compromise between the number of analysed stars and the errors 
on the trends. All the quality cuts on distance and atmospheric 
parameter errors remain similar to the ones of Sect. 2.3. 

4.4. Properties within 7.5 < R< 8.5 kpc 

We start by describing our results in the Solar suburbs (7.5 < 
R < 8.5 kpc) at two different height-bins from the Galactic 
plane: jZj < Ikpc and 1 < \Z\ < 2kpc. Figure 12 shows 
the matching of the recovered trends with the data-points at dif¬ 
ferent S/N cuts. Fig. 13 represents the [Mg/Fe] dispersions as 
a function of [Fe/H], and finally middle plots of Fig. 14 and 
Fig. 15 show the derived relative numbers of the two sequences 
close and far from the plane, and the mean required populations 
to fit the data. For both heights, the two populations are recov¬ 
ered nicely, with no particular indication of a merging of the se¬ 
quences up to [Fe/H] « -1-0.2. Indeed, the mean number of pop¬ 
ulations above [Fe/H] > 0 (second column of Fig. 15) is most 
of the time greater than or close to 1.5, however with a noisy and 
decreasing trend indicating the difficulty of fitting robustly our 
data. 

As mentioned in the previous Section, the data do not al¬ 
low us to recover reliably the position of a ‘knee’ in the [a /Fe] 
trend. However, we do recover the slopes of the decrease of the 
[Mg/Fe] — [Fe/H] for each of the low—a and high—a sequences. 
These slopes depend on the histories of star formation and gas 
accretion/outflow (e.g.: Snaith et al. 2014; Nidever et al. 2014, 
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Input 




[Fe/H] 



-1.5 -1.0 -0.5 0.0 0.5 

[Fe/Fl] with errors 



[Fe/H] 


Fig. 10. Top; adopted trends for the mock data used to test our fitting procedure (left) and effect on the mock data of introducing errors in 
[Fe/H] and [a/Fe] as the observations (right). Bottom: Output of our fitting procedure for the high—a (red crosses) and low—a (green squares) 
populations performed on the noisy mock data. Over-plotted in black are the error-free points used to create the mock catalogue. The bottom 
right plot shows the average of the number of populations needed to better fit the data based on the log-likelihood ratio tests of Eq. 6. When the 
red (green) line is below 1.5, this indicates that fitting the data with only a thick (thin) disc rather than two components simultaneously is more 
significant for more than half of the Monte-Carlo realisations. 


see also Sect. 5). The measured slopes are reported in Table 1 
and are discussed in the following subsections. 

4.4.1. Close to the plane 

For the stars close to the plane, in the metallicity interval 
[—0.9,0.0] dex we find (see Table 1): 

;xthi„ = (-0.20 ± 0.04) X [Fe/H] -f (0.00 ± 0.02) (7) 

/ithick = (-0.31 ± 0.05) X [Fe/H] -f (0.12 ± 0.02) (8) 


0<lZl<1 kpc, 7.5<R< 8.5 kpc, S/N> 25 



1<IZI<2 kpc, 7,5<R< 8-5 kpc, S/N> 25 



These slopes are robust to different S/N cuts (see Fig. 11; 
note we adopted S/N> 25), and are found to fit well the trends of 
Adibekyan et al. (2013) (grey points in Fig. 11, in the panel rep¬ 
resenting the Solar suburb). In addition, we found that increasing 
the S/N threshold emphasises the gap between the two chemi¬ 
cal sequences by depopulating the region between the chemical 
sequences (see Fig. 12). The prominence of this gap is there¬ 
fore compatible with the results of Bensby et al. (2014); Recio- 
Blanco et al. (2014), we note, however, that in order to establish 
whether this gap is true or just a trough, depopulated due to low 
number statistics, additional data will be needed (e.g. the future 
data-releases of Gaia-ESO). 

Our results show a steeper slope for the [Mg/Fe] sequence 
of the thick disc compared to the thin disc. As we will discuss 
in the following sections, a steeper slope implies higher ratio of 
past to present SFR (provided there is a fixed stellar IMF and 
Type la supernova delay distribution function for both discs). 
Above Solar metallicities (where the determination of the slope 
is not performed), we notice a flattening of the sequence for the 
thin disc (see Fig. 12) that could be compared to the “banana 


Fig. 13. [Mg/Fe] dispersions as a function of [Fe/H[, for the thick disc 
stars (in black) and the thin disc stars (in blue), close to the Galactic 
plane (left) and far from it (right). 


shape” quoted in Nidever et al. (2014, their Fig. 10), also seen in 
Bensby et al. (2014, their Fig. 15). 

We find that the a—enhancement dispersion, i.e. the thick¬ 
ness of the chemical sequence followed by either one of 
the discs, is of the order of ^ 0.08 dex and 

F’’[Mg/Fe]thick ~ 0.07 dex, in agreement with Mikolaitis et al. 
(2014), obtained using the Gaia-ESO iDRl. We note, however, 
that below [Fe/H] < —0.7 dex, the thin disc dispersion increases 
to 0.1 dex; this is due to low number statistics, as our test on the 
significance of the fit indicates that below this metallicity, one 
component model is on average a better fit to the data. 

For the considered local sample, the relative proportion of 
the low—a over high—a stars varies from 30 per cent to 80 
per cent in the iron abundance range between [—0.8, 0.0] dex, 
reaching 50 per cent at [Fe/H]~ —0.4 dex (middle plots of 
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Fig. 11. Trends of the [Mg/Fe] — [Fe/H] sequences for the thick disc (plain line) and thin disc (dashed line), for different signal-to-noise (S/N) 
threshold limits close to the plane (0 < \Z\ <1 kpc, left column) and far from it (1 < jZj < 2 kpc, right column). The top, middle and bottom 
rows show the Gaia-ESO iDR2 results for the inner Galaxy, Solar suburb and outer Galaxy, respectively. The over-plotted grey points in the panel 
7.5 < R < 8.5 kpc and 0 < \Z\ < Ikpc are the a—abundance measurements of Adibekyan et al. (2013), obtained for a very local sample 
(d < 150 pc). No significant differences in the derived trends are found as a function of S/N, which is indicative of the robustness of our fitting 
method. 


Table 1. Slopes of the [Mg/Fe] — [Fe/H] relations for thin and thick disc stars obtained for S/N> 25. 


\Z\{kpc) < 1 1 < |Z|(kpc) < 2 




a 

b 

a 

b 

Thin disc 

R[ kpc) < 7.5 
7.5< i7(kpc) < 8.5 
R{ kpc) > 8.5 

-0.14 ±0.07 
-0.20 ±0.04 
-0.23 ±0.03 

0.03 ±0.04 
0.00 ±0.02 
0.00 ±0.02 

-0.16 ±0.04 
-0.17 ±0.05 
-0.14 ±0.08 

0.06 ±0.02 
0.04 ±0.02 
0.03 ±0.04 

Thick disc 

7Z( kpc) < 7.5 
7.5< i7(kpc) < 8.5 
R{ kpc) > 8.5 

-0.27 ±0.05 
-0.31 ±0.05 
-0.28 ±0.07 

0.14 ±0.02 
0.12 ±0.02 
0.13 ±0.03 

-0.25 ±0.05 
-0.21 ±0.05 
-0.19 ±0.15 

0.19 ±0.02 
0.19 ±0.02 
0.14 ±0.07 


Notes. The a and b coefficients are defined for —0.9 < [Fe/H](dex) < 0, as: [a/Fe]^ a x [Fe/H]+6. 


Fig. 14). Below [Fe/H] ^ —0.5 dex, we find that the propor¬ 
tion of thin disc to thick disc stars is slowly decreasing down to 
[Fe/H] « —0.8 dex. This result is in agreement with the anal¬ 
ysis of Sect. 3.3, where we investigated the change in the rela¬ 
tive proportions of the different stellar populations by analysing 
the variations in 9V0/c)[Mg/Fe]. Indeed, recall that in Sect. 3.3 
we found that below —0.5 dex, and at least down to —0.8 dex, 
the value of dYd\Mg/¥e\ varied only slowly as a function 
of metallicity (albeit with large error bars) which is consistent 
with only a small variation in the relative number of stars of 
each disc within this metallicity range. This therefore also im¬ 
plies that the shape of the metallicity distribution functions of 
the discs are significantly non-Gaussian (see also Hayden et al. 
2015). We note that these results are robust to S/N selections 


and should not be dependent of the selection function, assuming 
there is no bias in [a/Fe] at fixed [Fe/Hj. For the most metal- 
poor bin ([Fe/H] = —0.9 dex), we find that the proportion of 
the low—a stars increases up to 45 per cent. Combined with the 
significance of the fit, and the increase of the [Mg/Fe] disper¬ 
sion (see above), this indicates that this increase is not real, and 
is rather due to low number statistics. 

Finally, we notice that both sequences are seen at all the 
metallicities, from —0.8 dex up to super-Solar values. However, 
the fact that the relative proportion trend reverses for [Fe/H] > 
0.1 (middle-top panel of Fig. 14) suggests that the number of 
stars in these bins is too small to perform a robust likelihood 
estimation and therefore separating the two populations when 
the sequences merge is challenging (as already mentioned in 
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Fig. 12. [Mg/Fe] as a function of [Fe/H] for the stars in the Solar neighbourhood, closer than Ikpc from the Galactic plane (left panels) and 
between 1 kpc and 2 kpc from the plane (right panels). Different selections, according to the signal-to-noise ratio of the spectra, are shown at 
the top (S/N> 10) and the bottom (S/N > 35). The colour-code corresponds to the azimuthal velocity of the stars. The dashed blue (black) line 
corresponds to the best fit model of the thick (thin) disc. The black ‘-I-’ sign is located at [Mg/Fe] = 0 and [Fe/H] = 0 dex. 


Sect. 4.3). This is also indicated from the plot of Fig. 15, where 
the significance of having only a thin disc rather than two popu¬ 
lations (blue line), is noisy above solar metallicities. 

We now assess to what extent the derived proportions are in 
agreement with the expected metallicity and spatial distribution 
functions of the discs. To do so, we model the expected number, 
N, of stars of a given disc as a function of Galactocentric radius, 


R, distance from the mid-plane, Z, and metallicity, [Fe/H], as: 


N{R,Z, [Fe/H]) 


/ ■ exp 


hz 


R — R{} 

hR 


•MDF, 


(9) 


where / is the local density of stars belonging to each of the 
discs, and (hz,hR) are the scale-height and scale-length. The 
MDF is the adopted metallicity distribution function, defined as: 
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0<lZl<1 kpc, 6.0<R< 7.5 kpc, S/N> 25 0<lZl<1 kpc, 7.5<R< 8.5 kpc, S/N> 25 0<lZl<1 kpc, 8.5<R<10.0 kpc, S/N> 25 



1<IZl<2 kpc, 6.0<R< 7.5 kpc, S/N> 25 1<IZl<2 kpc, 7.5<R< 8.5 kpc, S/N> 25 1<IZl<2 kpc, 8.5<R<10.0 kpc, S/N> 25 



Fig. 14. Relative proportion of the high—a (in black) and low—a (in blue) populations as a function of [Fe/H], for increasing Galactocentric radii 
(from left to right) and distances from the Galactic plane (upper panels: close to the plane, lower panels: far from the plane). Plain thick lines 
represent measurements derived from observations, smoothed with a boxcar average taking into account the closest neighbour, whereas thin red 
and grey lines are derived from the toy model described in Table 2, with the scale-heights and scale-lengths of Juric et al. (2008) and Bovy et al. 
(2012b), respectively. 
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Fig. 15. Number of required populations, averaged over 100 Monte-Carlo realisations, as a function of metallicity. For each realisation, the sig¬ 
nificance of the distribution is evaluated, and the likelihood of requiring one or two populations is evaluated. Blue (black) line represents the 
results for the low—a (high—a) population. A value larger than 1.5 indicates that the low—a (high—a) population cannot describe by itself the 
observed trends. 


3 

MDF = ^ ai ■ exp 

i=l 

It is obtained as the sum of three Gaussians, associated with 
the metal-poor, intermediate metallicity, and metal-rich regimes, 
in order to reproduce skewed distributions. The factor is the 
weight assigned to each of these metallicity regimes, while Mi 
and CTi are the mean metallicity and dispersion of the associated 
Gaussian. We assume that the MDF has a fixed shape, but can 
shift to lower or higher mean metallicities as a function of R and 
Z. We therefore define Mi, as follows: 

M, = fj,, + Z- a[Fe/H]/5Z + (R - Rq) • (9[Fe/H]/(9R, (11) 


where R is expressed in kiloparsecs and fii is the mean of the 
metal-poor, intermediate and metal-rich regimes, as measured at 
the Solar neighbourhood. The input values for our simple expec¬ 
tation models are summarised in Tables 2 and 3. 

The adapted parameters are not chosen in order to fit the 
Gaia-ESO data, but are a compilation derived from the literature, 
based on discs defined either by kinematics, chemical abundance 
or star-counts. Therefore, the parameters of Table 2 roughly re¬ 
produce the skewness of the MDFs for the metal-weak regime 
of the thin and thick discs (Wyse & Gilmore 1995; Kordopatis 
et al. 2013b) and the super-solar metallicity stellar distribution 
of RAVE (Kordopatis et al. 2015). The vertical and radial metal¬ 
licity gradients of the thin disc are adopted from Mikolaitis et al. 
(2014) and Gazzano et al. (2013), respectively. Finally we adopt 


([Fe/H]-M,)^ 




( 10 ) 
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Table 2. Input values for the metallicity distribution functions of the 
discs. 

thin disc thick disc 
7 085 0T4 

ai(dex) 0.90 0.80 

/ri(dex) 0.0 —0.5 

cri(dex) 0.08 0.35 

a2(dex) 0.1 0.15 

/i2(dex) —0.3 —0.8 

cr2(dex) 0.3 0.5 

a3(dex) 0.02 0.05 

/r3(dex) +0.2 +0.2 

0-3 (dex) 0.5 0.5 

9[Fe/H]/9R(dexkpc-i) -0.05 0.00 

a[Fe/H]/9Z (dexkpc-i) -0.06 0.00 


Table 3. Adopted scale-lengths and scale-heights of the discs. 



Thin disc 

Thick disc 


hR 

hz 

hR 

hz 


(kpc) 

(kpc) 

(kpc) 

(kpc) 

Juric et al. (2008) 

2.6 

0.30 

3.6 

0.9 

Bovy et al. (2012b) 

3.8 

0.27 

2.0 

0.7 



the scale-heights and scale-lengths of June et al. (2008, defined 
morphologically) and Bovy et al. (2012b, dehned chemically, 
see their Tables 1 and 2), to illustrate how the star counts are 
expected to change when considering a thick disc that is more 
extended or less extended than the thin disc. 

The expected proportions of thin and thick discs are plotted 
in Fig. 14. One can see that for the sample in the range 7.5 < 
R < 8.5 kpc, close to the plane (middle top panel), the expected 
trends are not matching correctly the observed ones: the model 
of Juric et al. (2008) over- (under-) estimates the thick (thin) 
disc relative proportions, whereas the inverse is noticed for the 
Bovy et al. (2012b) model. We note, however, that modifying 
the Bovy et al. (2012b) model, with a thin disc slightly more 
compact {hfi = 3.6 kpc, hz = 0.25 kpc) and a thick disc with 
a larger scale-height {hz = 0.9 kpc), gives a trend in a better 
agreement with the observations. 

4.4.2. Far from the plane 

For the stars being between 1 and 2 kpc from the Galactic plane 
(bottom plots of Figs. 12 to 15), we hnd that the slopes of the 
mean a—abundances are slightly flatter, compared to the derived 
values closer to the plane (/ithin = (—0.17 ± 0.05) x [Fe/H] + 
0.04, /ithick = (-0.21 ± 0.05) X [Fe/H] + 0.19). In addition, 
the mean a—enhancement of the discs at a given metallicity is 
increased by ^ 0.04 dex for the thin disc and ~ 0.07 dex for 
the thick disc (Fig. 17), in agreement with the vertical [a/Fe] 
gradients in the discs found by Mikolaitis et al. (2014). 

In addition, despite the fact that the selection function of 
Gaia-ESO is not taken into account (see Sect. 2), we hnd that 
the relative weight of the two populations changes, with the rel¬ 
ative importance of the thick disc being higher, as expected (e.g.: 
Minchev et al. 2014). At Solar metallicities, at these heights 
above the plane, the thick disc represents 50 per cent of the 
stars. Compared to our star-count models, the agreement with 
the Bovy et al. (2012b) model is very good, in particular for 
the metal-poor part ([Fe/H] < —0.4), however with again an 


Fig. 16. Trends at different Galactocentric radii for given bins of dis¬ 
tances from the Galactic plane. Plain and dashed lines are associated to 
the high—a and low —q populations, respectively. 


over-estimation of the thin disc relative proportion at the higher 
metallicities. 

Finally, the [a/Fe] dispersions show mild trends (Fig. 13): 
the thick disc’s one decreases from 0.07 dex at low metallicities 
to 0.05 dex at super-solar values, as does the thin disc’s, with 
however an important increase for metallicities below —0.6 dex. 
Nevertheless, this increase is only due to the low number of stars 
available to fit for the low—a sequence, as shown by the relative 
proportion in Fig. 14 and the signihcance of the ht for metallici¬ 
ties below [Fe/H] ~ —0.6 dex (Fig. 15). 

4.5. Investigation beyond 7.5 < R< 8.5 kpc 

The same htting procedure as for the Solar suburbs has been 
repeated for radii intervals between 6.0 < R < 7.5 kpc and 
8.5 < R < 10.0 kpc. The bin sizes are purposely wider than 
the ones at the Solar suburb, in order to have enough stars to 
analyse, and in order to minimise the impact of having an inho¬ 
mogeneous sky coverage (see Fig. 3). Figure 14 shows how the 
relative proportion of each population changes as a function of 
the region in the Galaxy. One can see that the proportion of the 
thick disc decreases with R, both close and far from the plane, 
suggesting that the thick disc is more radially concentrated than 
the thin disc, in agreement with previous surveys (e.g.: Bensby 
et al. 2011; Cheng et al. 2012; Bovy et al. 2012b; Anders et al. 
2014; Recio-Blanco et al. 2014; Mikolaitis et al. 2014), and that 
the metal-weak tail of the low—a population does not extend be¬ 
low —0.6 in the 6.5 < R < 7.5 kpc ring. Our star-count models 
achieve a fitting of the relative proportions in different ways: for 
the range 6.0 < R < 7.5 kpc, the model of Juric et al. (2008) 
(thick disc longer than the thin disc) seems more appropriate, 
however, the Bovy et al. (2012b) model, which is assuming a 
shorter scale-length for the thick disc, is in better agreement for 
the stars located at the largest radii. This disagreement could in- 
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Fig. 17. Trends at different distances from the Galactic plane: 6 < R < 
7.5 kpc (top), 7.5 < R < 8.5 kpc (middle) and 8.5 < R < lOkpc 
(lower panel). 


dicate a varying scale-length or scale-height with R for either or 
both of the discs, as recently suggested by Minchev et al. (2015). 

The [Mg/Fe] — [Fe/H] slopes as a function of [Fe/H] at the 
three Galactic regions are shown in Fig. 16, and summarised in 
Table 1. Overall we find that as a function of R, the thick disc has 
similar trends, whereas the thin disc slopes become marginally 
steeper (from —0.14 ± 0.07 at i? < 7.5 kpc to —0.23 ± 0.03 
at i? > 8.5 kpc). In addition the slopes become, on average, 
flatter as a function of Z, by 0.03 ± 0.05 for the thin disc and by 
0.07 ± 0.04 for the thick disc. 

The result from the Gaia-ESO survey obtained here, that the 
high-cr thick disc sequence varies little spatially, is consistent 
with the conclusions from the APOGEE data for red clump stars 
published recently by Nidever et al. (2014). Those authors found 
that the overall spatial variations of the [a /Fe] — [Fe/H] sequence 
of the high—a stars do not exceed 10 per cent. 

Finally, we find that the mean a—enhancement for both discs 
are increased higher from the plane, at all Galactocentric radii. 
This is indicative of a vertical gradient in a—elements in the thin 
disc, in agreement with the value of 4-0.04 dexkpc”^ measured 
in Mikolaitis et al. (2014). 


5. Discussion: Impiications for how the discs 
evoived 

The steepness of the [a/Fe] — [Fe/H] slopes contains informa¬ 
tion about how discs have evolved, encoding mainly informa¬ 
tion on the time evolution of their Star Formation Rates (SFR) 
and gas in/out flows, assumed to regulate the SFR. Typically, the 
[a/Fe] decline with increasing [Fe/H] is caused by the time de¬ 
lay distribution that the first supemovae of type la (mostly iron 
producers) have after the core-collapse supernovae (mainly pro¬ 
ducers of a—elements). After the adoption of an IMF, binary 
fraction and delay time distribution, the relative present to past 
SFR can be translated into a relative rate of core-collapse to Type 


la supernovae. In turn, this defines a slope of [a/Fe] — [Fe/H] 
for the ISM and hence for newly formed stars, modulo possible 
elemental-dependent gas in/out flows. 

In the case that most of the stars are formed during a star- 
burst event, followed by a negligible star formation, then the 
slope of the decline - tracing the abundance of the gas - is max¬ 
imal: the [a/Fe] — [Fe/H] slope connects the last star formed 
before the end of the starburst and the hrst one to form after the 
gas has been enriched by SNIa ejecta. On the other hand, flatter 
slopes, i.e. slower declines, imply more important contributions 
from massive stars through core-collapse supemovae and thus 
more extended star formation. A typical slope for a closed-box 
star-forming system with a slowly decreasing SFR as inferred 
for the solar neighbourhood and standard Type la delay time 
distribution is around —0.3 (see e.g. Fig. 1 of Gilmore & Wyse 
1991). The flattening of the slope as the star-formation efficiency 
(defined as the reciprocal of the duration of star formation) de¬ 
creases is illustrated in the top panel of Fig. 15 of Nidever et al. 
(2014). The middle panel of that figure in Nidever et al. (2014) 
illustrates how changing the rate of gas outflow can also change 
the slope. 

Our work finds steeper slopes for the thick disc (« —0.28 x 
[Fe/H]) compared to the thin disc, with negligible spatial varia¬ 
tions as a function of R and Z. On the other hand, the thin disc 
slopes show a variation as a function of R, with the flatter slope 
being in the inner Galaxy (—0.14 x [Fe/H]) and the steeper in 
the outer Galaxy (—0.23 x [Fe/H]). 

Edvardsson et al. (1993) undertook a pioneering study of 
kinematics and elemental abundances, for a sample of 189 F and 
G dwarfs observed at the solar neighbourhood. They investigated 
the manner in which the [a/Fe] abundance of a star in a given bin 
of [Fe/H] depends on the star’s mean orbital radius, using this 
as a proxy for birth radius. These authors found that stars from 
the outer Galaxy (mean orbital radius R^ > 9 kpc) follow an 
[a/Fe] — [Fe/H] relation that falls slightly below that followed by 
stars with mean orbital radius close to the Solar Circle. Further, 
they found that for stars from the inner Galaxy {Rm < 7 kpc), 
the [a/Fe] values split into two groups; one at [a/Fe] = 0.2 dex 
for [Fe/H] < —0.4 dex, and one at [a/Fe] = 0.03 dex for 
[Fe/H] > —0.4 dex. The conclusion of their analysis was that 
the variation of the [a/Fe] — [Fe/H] sequences is evidence of 
inside-out formation of the disc, as the inner disc region had a 
higher past to present SFR, and the outer parts form more slowly. 

On the other hand, the two infall model of Chiappini et al. 
(1997), that also forms the thin disc in an inside-out fashion 
however with contribution of accreted pristine extra-galactic gas 
once the thick disc is formed, predicts opposite trends: the inner 
Galaxy has flatter [a/Fe] — [Fe/H] trends, as a consequence of 
more important quantity of infalling gas and therefore more ex¬ 
tended star formation history. Nevertheless, the difference in the 
trends that is predicted by this model is rather small and is only 
noticeable over a large radial range (from 4 kpc up to 16 kpc, see 
Chiappini et al. 1997, their Fig. 14). It is clear that the results are 
quite sensitive to the details of the adopted model. 

Our results, based upon an analysis using the observed po¬ 
sitions of the stars'^ And that the slopes of the [a/Fe] — [Fe/H] 
sequences are flatter for thin-disc stars in the inner Galaxy, com¬ 
pared to the outer parts. However, the radial range that is spanned 
by the Gaia-ESO targets, combined with the fact that the spatial 


As opposed to an estimate of the mean orbital radius used 
by Edvardsson et al. (1993). The radial excursions associated with 
epicyclic motions for thin disc stars at the Solar neighbourhood with 
Gv ~ 30 km are of the order 1 kpc (Binney & Tremaine 2008). 
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variations are below 2a, cannot allow us to discard either of the 
models. 

On the other hand, we find almost no variations (differences 
are below Icr) of the trends for the thick disc, in agreement with 
recent results from Nidever et al. (2014) using APOGEE data. 
To explain the invariance of the thick disc sequence, Nidever 
et al. (2014) suggested the existence of a large-scale mixing in 
the early “turbulent and compact” disc. We note, however, that 
in such a disc, the velocity dispersion of the gas should be high, 
in order to mix efficiently the metals from the inner Galaxy with 
regions at i? ss 5 kpc. This could potentially be incompatible 
with the observed velocity dispersion of the thick disc stars that 
formed from that gas, unless they were formed during dissipa- 
tional collapse. Such a collapse would imply vertical gradients 
in metallicity, a—elements and azimuthal velocity, all being still 
a matter of debate in the literature (see Mikolaitis et al. 2014, for 
an overview of the published values), though being compatible 
with our data. More theoretical work is needed to quantify the 
expected amplitude of any vertical gradient. 

Indeed, we found that the a—enhancement of the popula¬ 
tions increase as a function of the distance from the mid-plane. 
Eor the thin disc, this vertical gradient is a signature of the age- 
velocity relation; older stars (the more a—enhanced) have larger 
random velocities than newly born stars, plausibly due to inter¬ 
actions with gravitational perturbations such as spiral arms and 
Giant Molecular Clouds. On the other hand, the saturation of the 
thin disc age-velocity dispersion at cr^ ~ 20 km (Nordstrom 
et al. 2004; Seabroke & Gilmore 2007), well below the disper¬ 
sion of the thick disc, rules out such an explanation for the verti¬ 
cal gradient in the thick disc. In addition, such a vertical gradient 
is difficult to explain by a single accretion event that formed the 
thick disc from the accreted stars, as for example suggested by 
Abadi et al. (2003). A gradient as observed would, however, nat¬ 
urally arise in the scenarios where the thick disc is formed (i) 
from the stars formed during dissipational settling (e.g. Brook 
et al. 2004; Minchev et al. 2013; Bird et al. 2013) (ii) by the 
dynamical heating of a pre-existent thin disc (with pre-existent 
vertical [a/Fe] gradient resulting from an age-velocity disper¬ 
sion relation) due to minor mergers (e.g.: Wyse & Gilmore 1995; 
Villalobos & Helmi 2008; Minchev et al. 2014), or (Hi) by a thick 
disc formed by radial migration (Schonrich & Binney 2009). 
We note, however, that the simulations of thick disc formation 
through stellar radial migration (e.g.; Loebman et al. 2011) have 
difficulties reproducing many of the observed thick disc proper¬ 
ties (Minchev et al. 2012). 

5 . 1. The low-metallicity extent of the thin disc 

The extent of the thin disc sequence to low metallicities and its 
overlap in iron abundance with that of the thick disc provides an 
important constraint on the evolutionary sequence of the discs. 
The low-metallicity limit is of particular interest in models in 
which there is a hiatus in star formation between thick and thin 
discs, with inflow of metal-poor gas, and good mixing, to re¬ 
duce the mean enrichment of the proto-thin-disc material (e.g. 
Haywood et al. 2013; Snaith et al. 2014, see also discussion 
in Wyse & Gilmore 1995). Obviously the lower the metal-poor 
limit of the thin disc, the larger the mass of metal-poor gas that 
must be accreted. Our work has shown that, close to the plane, 
the thin disc extends at least down to [Fe/H] ~ —O.Sdex in 
the Galactocentric radial range 7.5 < R < 10 kpc and at least 
down to [Fe/H] ~ —O.Odex in the range 6.5 < R < 7.5 kpc. 
A lower limit of —0.6dex also seems to be reached far from 
the plane (1 < \Z\ < 2 kpc), at all the investigated radii. The 


presence of thin-disc stars below these metallicities is less se¬ 
cure. Indeed, given the uncertainties on the relative proportion 
of the populations the uncertainties on the [a/Fe] dispersions as 
well as the uncertainties on the results of the log-likehood ratio 
tests, we cannot state with certainty that the thin disc, as defined 
by the low—a elemental abundance sequence, extends down to 
[Fe/H] ~ — 1 dex. A more robust determination of the propor¬ 
tion of thick- to thin-disc stars requires that the selection func¬ 
tion of the Gaia-ESO Survey be corrected for, and this is beyond 
the scope of the present paper. We note that when the thick and 
thin discs are defined through their vertical kinematics or ver¬ 
tical extent from the mid-plane, very similar results to those of 
this paper are obtained for the metal-poor extent of the thin disc 
- down to ~ —0.8 dex in iron abundance - and for the relative 
proportions of thick- and thin-disc stars below —0.4 dex (Wyse 
& Gilmore 1995). 

The metal-poor low—a stars could have formed after sig¬ 
nificant metal-poor gas infall into the disc, or, assuming inef¬ 
ficient mixing, in a relatively metal-poor region. Alternatively 
these stars could have formed in the outer Galaxy, where the gas 
was less-enriched, and migrated radially inwards. In the absence 
of good estimates of the individual stellar ages, it is difficult to 
distinguish these possibilities. 

We nevertheless searched within our Solar suburb sample for 
contamination by stars formed at other radii. This contamination 
could include stars indistinguishable by their kinematics (i.e. ra¬ 
dial migrators from corotation resonances with the spiral arms, 
see Sellwood & Binney 2002) and stars on apo- and pericenters 
currently visiting the solar region (which could also be radial 
migrators that changed their angular momentum at the Lindblad 
resonances with the spiral arms and thus are on more elliptical 
orbits, see Sellwood & Binney 2002). Since velocity can give us 
a partial hint on the origin of the stars, we used the slopes de¬ 
fined in Table 1, to select stars that have a [Mg/Fe] ratio below 
these trends for the thin disc, and above these trends for the thick 
disc. We then investigated how their azimuthal velocity distribu¬ 
tions compared. The stars residing between these cuts were then 
labelled as “intermediate—a” stars (i.e. potentially being inside 
the gap). 

The histograms of Fig. 18, obtained for the sample having 
S/N > 10, show that the low—a population has a broader ve¬ 
locity distribution at the metal-poor regime than at the metal-rich 
regime, for all the three investigated Galactic volumes. This is 
indicating that most (but not all) of the these metal-poor low—a 
stars have been kinematically heated by interactions at Lindblad 
resonances and that they are very likely close to their pericentre 
(i.e. their guiding centre is at larger radius than their observed 
one). Indeed, assuming a negative metallicity gradient for the 
interstellar medium (e.g. Genovali et al. 2014), it follows that 
many of these low metallicity stars would have been formed in 
the outer Galaxy. Furthermore, noticing that the velocity distri¬ 
bution of the low metallicity low—a stars is also shifted towards 
higher V^, this again implies that these stars are at their peri- 
galacticon, increasing their velocity as their Galactocentric ra¬ 
dius decreases. We note that on the other hand, the high—a pop¬ 
ulation exhibits the typical correlation between the azimuthal ve¬ 
locity and the metallicity of the thick disc (measured for a thick 
disc defined by metallicity only, e.g., Kordopatis et al. 2011, 
2013b,c). From that, we conclude that it does not show any par¬ 
ticularly odd kinematic behaviour for its low metallicity stars. 

Finally, in order to assess whether there is a true gap between 
the high—a and low—a chemical paths, we have computed the 
Kolmogorov-Smirnov (KS) probabilities, between the velocity 
distribution of the intermediate—a population and the ones of 
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- 1 .0<[Fe/H]<-0.4 ; R<7.5 kpc 



-0.4<[Fe/H]< 0.0 ; R<7.5 kpc 



-1.0<[Fe/H]<-0.4 ; 7.5<R<8.5 kpc 


-1.0<[Fe/H]<-0.4 ; R>8.5 kpc 



-0.4<[Fe/H]< 0.0 ; 7.5<R<8.5 kpc 




Fig. 18. V<^—histograms of the high —q (orange), low—a (red) and intermediate—a (green) stars at two different metallicity regimes: metal-poor 
(top) and intermediate (bottom) metallicity. The stars are selected to be close to the Galactic plane (\Z\ < 1 kpc) at three Galactocentric radial 
ranges (left panels: R < 7.5 kpc, middle panels: 7.5 < R < 8.5 kpc, right panels: R > 8.5 kpc). The Kolmogorov-Smimov probabilities that 
the intermediate—a histograms are drawn from the distributions as the low— or high—a histograms are shown in Table 4. 


Table 4. Kolmogorov-Smirnov probabilities between the azimuthal ve¬ 
locity distribution of the intermediate-a population and the ones of the 
low— and high—a populations. 
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[M/H] > -0.4 

3.2- lO-'^ 

1.0- lo-'^ 

R > 7.5 kpc 

[M/H] < -0.4 

0.42 

0.30 


[M/H] > -0.4 

0.51 

0.02 


the thin and thick discs (green, red and orange histograms in 
Fig. 18, respectively). These probabilities quantify the signifi¬ 
cance of two datasets being drawn from the same distributions. 
The results shown in Table 4 indicate that the intermediate—a 
velocity distribution is always more resembling to the low—a 
velocity distribution than the high—a one, therefore suggest¬ 
ing that there is a true gap between the sequences, rather than 
a trough. 

5.2. The high-metallicity extent of the thick disc 

Similar to the considerations of the low metallicity tail of the 
thin disc, the high-metallicity extent of the thick disc holds in¬ 


formation about the degree of interdependence in the evolution 
paths of the discs. Indeed, Haywood et al. (2013) and Bensby 
et al. (2014) suggested that at least part of the metal-poor end of 
the thin disc could have started forming from either accreted gas 
at the end of the star formation epoch of the thick disc and/or 
from gas expelled from the thick disc (see also Gilmore & Wyse 
1986, for an example of the latter type of model, or Chiappini 
et al. 1997 for the two infall model). 

We have demonstrated the presence of high—a stars up to 
Solar metallicities, and probably up to +0.2 dex. This result is 
in agreement with previous studies, (e.g.: Adibekyan et al. 2013; 
Bensby et al. 2014; Nidever et al. 2014, most recently); see also 
Bensby et al. (2007); Kordopatis et al. (2013b, 2015) for fur¬ 
ther evidence of metal-rich thick disc stars, dehned kinemati¬ 
cally. We hnd that the dispersion in [Mg/Fe] at hxed iron abun¬ 
dance for thick-disc stars is marginally smaller than the values 
derived for the thin disc. Assuming that this dispersion reflects 
initial conditions, the smaller dispersion is indicative of a better- 
mixed ISM for the thick disc compared to the thin disc. Possible 
reasons for that could either be accreted gas that would change 
the ISM composition, or stellar radial migration that would mix 
stars formed in different regions of the Galaxy. However, even if 
our results seem to indicate a separation between the sequences 
over the entire [Fe/H] range, this might not be the case for the 
most metal-rich end, since our procedure always tries to ht two 
components to the data. 
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Finally, we note that Martig et al. (2014) and Chiappini 
et al. (2015) have found by combining APOGEE abundances and 
astero-seismic ages, high—ct stars with ages younger than 6 Gyr, 
including in the super-Solar metallicity regime. Such young ages 
would imply that these stars were formed within the Galaxy, 
reflecting the complex star formation history and ISM mixing 
of the Milky Way in the inner most regions. The quality of our 
data does not allow us to derive reliable ages by isochrone fitting 
for our stars. We therefore cannot exclude the existence of such 
high—Of metal-rich stars in our sample that would not belong to 
the thick disc but rather in the thin disc. 

6. Summary 

We have investigated the stellar abundances and kinematics 
coming from the second internal data release of the Gaia-ESO 
Survey in order to characterise, in a statistical way, the high—a 
and low—a disc populations, commonly identified with the thin 
and thick discs. 

By studying the way the azimuthal velocity correlates with 
the a—abundances at a given metallicity, we evaluated the extent 
of the tails of the metallicity distribution function of the discs, 
without making any assumption on their shape. We found that 
the thick disc extends up to super-Solar metallicities ([Fe/H] « 
-1-0.2) and the thin disc down to [Fe/H] « —0.8. The kine¬ 
matic data shows, in agreement with a decomposition of the 
discs according to their a—abundances, that at the extended 
Solar neighbourhood, there is an almost constant relative frac¬ 
tion of thin disc and thick disc stars over the metallicity range 
-0.8 < [Fe/H] < -0.5. 

In a complementary approach, we also investigated the dis¬ 
tribution of stars in chemical space for different locations in the 
Galaxy. We found a mild steepening (with a significance below 
2a) of the decline of the [a/Fe] — [Fe/H] trend for the thin disc 
towards the larger Galactic radii. This could potentially suggest 
either a more extended star formation history in the inner Galaxy, 
or strong effects from radial migration, flattening the slope of the 
[a/Fe] — [Fe/H] relation over time. However, the relatively short 
radial range (6 < i? < 10 kpc) covered with the Gaia-ESO tar¬ 
gets does not allow us to put strong constraints or discard any 
model. Eor the thick disc, no variations are detected across the 
Galaxy, suggesting that the processes that formed or mixed the 
thick disc stars were very efficient. 

We found that the way the relative proportion of the disc 
stars varies across the Galaxy confirms the recently suggested 
hypotheses that the thick disc defined by a—abundance has a 
shorter scale-length than the thin disc. However, when compar¬ 
ing with toy models of thick and thin discs relative fractions hav¬ 
ing shorter or longer scale-lengths, we find evidence of varying 
scale-lengths and/or scale-heights with R. In a recent analysis, 
Minchev et al. (2015) suggested that in inside-out disc formation 
models, the thick disc (defined morphologically) is composed of 
the nested flares of mono-age populations. It was shown that flar¬ 
ing increases for older stellar samples, due to the conservation of 
vertical action when stars undergo radial migration. Since older 
populations have generally lower metallicities, such a flaring 
thin disc would populate mostly the regions at [Fe/H] < 0.4 dex 
above the plane (Fig. 14, bottom), thus decreasing the fraction of 
the chemically defined thick disc. 

Finally, by selecting the stars residing between the two 
chemical sequences, we have investigated their kinematics and 
found that they had typical thin disc velocities, suggesting that 
the thin and thick discs are following chemical paths separated 
by a true gap. However, the quality of our data does not allow us 


to establish firmly whether the two chemical sequences merge 
at super-solar metallicities. In future work, we plan to investi¬ 
gate part of the Gaia-ESO UVES sample (Smiljanic et al. 2014), 
with stellar parameters obtained from spectra at a resolution of 
R ^ 47 000, in order to investigate more robustly the metal-rich 
end of the discs. 
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